Renormalization group approach to the critical behavior of the forest fire model 
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We introduce a Renormalization scheme for the one and two dimensional Forest-Fire models in 
order to characterize the nature of the critical state and its scale invariant dynamics. We show the 
existence of a relevant scaling field associated with a repulsive fixed point. This model is therefore 
critical in the usual sense because the control parameter has to be tuned to its critical value in 
order to get criticality. It turns out that this is not just the condition for a time scale separation. 
The critical exponents are computed analytically and we obtain v = 1.0, r = 1.0 and v — 0.65, 
r = 1.16 respectively for the one and two dimensional case, in very good agreement with numerical 
simulations. 

PACS numbers: 64.60.Ak, 02.50.-r, 05.40.+j 



The Forest-Fire model (FFM) has been intro- 

duced as a possible realization of Self Organized Criti- 
cality (SOC) ||. The model formulated first by Bak and 
coworkers |lj] contains a tree growth probability p and 
fire spreading to nearest neighbours. In the limit of very 
slow tree growth, this model becomes more and more de- 
terministic [0 and shows spiral-shaped fire fronts. For 
this reason, the original model was modified by Drossel 
and Schwabl Q by introducing an ignition parameter /, 
the lightning parameter. This parameter is the proba- 
bility that during one time step a tree without burning 
nearest neighbours becomes a burning tree. In this case 
a critical behaviour, in the sense of anomalous scaling 

laws, is observed in the double limit 6 = - —> and 

... p . 

p — > 0. These two limits describe a double separation of 

time scales: trees grow fast compared to the occurrence 
of lightning in the system and forest clusters burn down 
much faster than trees grow. This means that each fire 
triggered by individual lightning does not overlap with 
other fires, thus clusters destroyed by fire are well de- 
fined objects. The critical state is characterized by a 
power law distribution P(s) ~ s~ r of the forest clus- 
ters of s sites and the average cluster radius scale as 
R ~ 6~ v '. Since this critical state is reached indepen- 
dently on the initial conditions and for a wide range of 
parameter values, it is called self-organized critical state. 
This statement however is ambiguous in that Self Orga- 
nized Criticality refers to the tendency of large dynamical 
systems to evolve spontaneously in a critical state with- 
out the fine tuning of any critical parameter. From this 
point of view, the role played by the parameter 6 is not 
clear. In fact it seems an effective relevant parameter, as 
the reduced temperature in thermal phase transitions, in 
that it allows criticality just for its critical value = 0. 

Usually the limit — > is referred to the existence of 
a slow driving of the system or, in other words, a time 



scale separation, also present in the definition of SOC 
models. It is worth to stress that in the SOC models the 
requirement of a slow driving of the system is essentially 
different. In this case, in fact, the time scale separation 
just affects the lower cut-off of the avalanche size dis- 
tribution, the upper cut-off being affected by finite-size 
effects ||. In the FFM the parameter 6, on the contrary, 
affects directly the upper cut-off Q and then it seems to 
play, in the language of ordinary critical phenomena, the 
role of a relevant parameter. 

In the context of Renormalization Group (RG) meth- 
ods, a self-organized critical phenomenon can be viewed 
as a critical phenomenon in which all the parameters 
are irrelevant, namely the scale invariant dynamics cor- 
responds to an attractive fixed point in the parameter 
space (phase space). To address the study of intrinsi- 
cally critical phenomena, we have developed the Fixed 
Scale Transformation approach to fractal growth and re- 
cently we have introduced an RG scheme of novel type 
to study sandpile models [§,1). 

In this letter we follow the same ideas in order to study 
the FFM including the ignition parameter /. This allows 
us to clarify the role of the critical parameters and p 
and the nature of the critical state. In fact, we obtain 
that is the relevant control parameter, so that FFM 
belongs more to the rich domain of ordinary critical phe- 
nomena than to self-organized criticality. In addition we 
are able to compute analytically the critical exponents 
characterizing the model in d — 1 and d = 2. 

For the sake of clarity we recall briefly the dynami- 
cal rules of the model we are going to discuss. At each 
time step the lattice is updated as follows: i) a burning 
tree becomes an empty site; ii) a green tree becomes a 
burning tree if at least one of its neighbors is burning; 

iii) at an empty site a tree can grow with probability p; 

iv) a tree without burning nearest neighbors becomes a 
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burning tree with probability /. 

Starting with arbitrary initial conditions, the system 
approaches, after a short transient, a steady state the 
properties of which depend on the parameter values. Let 
Po, pi and pi be the mean densities of empty sites, of trees 
and of burning trees in the stationary state. These prob- 
abilities describe in full generality the stationary proper- 
ties of the model, while, as for the dynamical properties, 
we will refer to a phase space defined by the parameters 
p and 8 = £ . In such a scheme p and 6 parameterize the 
dynamics of the system. 

We should now extend the characterization of the sta- 
tionary and dynamical properties of the system at a 
generic scale — bo2 k by considering coarse-grained 
variables pW, and p^ = (p ( k) , pf\ p ( 2 k) ). In the 
following we will mainly refer to the two-dimensional 
case but all the procedure can be drawn also in the one- 
dimensional case. The stationary properties of coarse 
grained variables are defined as follows (Fig.l): A cell of 
size b 1 -^ is considered as "green" if it is spanned from left 
to right by a connected path of green sites at scale &C C ~ 1 ). 
On the opposite a cell is empty if it is not spanned by 
a connected path of green sites. Finally, we consider a 
cell as burning if it contains at least one burning tree. 
In the latter case the spanning condition is not necessary 
because the fire spreads automatically to nearest neigh- 
bors sites. Consequently, the dynamical properties of the 
coarse grained description of the system will be charac- 
terized by the variables p^ and 6^ = The first 
variable express the probability that an empty cell of size 
&( fe ) becomes a green tree. The second one is defined using 
the probability that a green cell of size b^ becomes 
a burning cell because of the arrival of a lightning event. 

We now proceed to define a renormalization transfor- 
mation for the dynamical variables. We will use a cell-to- 
site transformation on the square lattice, in which each 
cell at scale b is formed by four sub-cells at scale 6/2. 
Every cell at the larger scale is then characterized by a 
number of green and empty sub-cells ranging from one to 
four. The relative weight of each green configuration is 
given by the probability to have the corresponding num- 
ber of green and empty sub-cells 

W a {p)=n aP «p\- a 1 (1) 

where a is the number of empty sub-cells and n a is a nor- 
malization factor that takes into account the multiplicity 
of the configurations. In order to define a renormaliza- 
tion transformation we start with an empty or green cell 
configuration at scale b and study how it evolves using 
the dynamical rules of the model. Here we consider a 
transformation defined by the following rules: i) Every 
series of tree growth processes at scale 6/2 that span an 
empty cell at scale 6 in the horizontal direction is renor- 
malized in the growth probability p at scale 6. ii) Every 
lightning process at scale 6/2 that affects a green cell of 
scale 6 contributes to the renormalization of the light- 
ning probability /. The spanning rule implies that only 



tree growth processes extending over the size of the new 
length scale contribute to the renormalized dynamics. 
Moreover, it ensures the connectivity properties of the 
green sites in the renormalization procedure. An exam- 
ple of such a renormalization procedure for the lightning 
probability is shown in fig. 2 for the simplest calculation 
scheme; i.e. 2x2 cells with left-right spanning condition.. 
We have two possible processes depending on whether 
both or just one green sub-cell is hit by a lightning event. 
The first case occurs with probability 2/( fe )(l - /(*>), the 
latter with probability (/ ( - fc - ) ) 2 . By summing these proba- 
bilities one obtains, the renormalization equation for the 
configuration with a = 2. We can write a renormaliza- 
tion equation for each configuration a corresponding to a 
green site, and, averaging over all the configurations we 
obtain: 

f^ = f2W a fi k+1) (2) 

a=2 

with 

f { 2 k+1) = / (fe) ( 2 - / (fe) ) 

/ 3 (fe+1) = /«(3-3/W + (/ w ) 2 ) (3) 
M+i) = j(fc)(4_gj(fc) + 4(y(fc))2 _ (/(*0)3) 

In an analogous way we can write the renormalization 
equations for p( fe + x ). Let us consider an empty cell at 
scale 6 of type a = 2 (see Fig. 3). After an updating step, 
the cell will become green if at least one of the two empty 
sub-cells has become green. This occurs with probability 
(p( k )) 2 and 2p( fc )(l — p^), corresponding to the fact that 
both or only one sub-cell become green, respectively. In 
a similar way one can also write a renormalization equa- 
tion for each configuration a corresponding to an empty 
site, and finally averaging over all the configurations we 
obtain 

p( k + 1 ) = j2W*p£ +1) (4) 

with 

Po +1) = (p (fe) ) 2 (2-(pW) 2 ) 

p (k+i) =p (fc) (!+„(*) -( P W) 2 ) (5) 

From the previous equations it is then straightfor- 
ward to obtain the renormalization equation = 
j{k+\) jpik+i) _ These renormalization equations are not 
yet complete because the statistical weights W a are func- 
tion of the the density vector p k . In fact, in order to de- 
scribe the stationary critical state it is necessary to couple 
the dynamics to a stationarity condition that gives the 
renormalization equations for the density vector. This 
scheme is similar to that used in ||. The stationarity 
condition is obtained from the master equations for the 
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density vector in the mean field regime by imposing the 
asymptotic equilibrium condition (t — ► oo) [p|,^o|, 



(l-p( fe) )a (fc) /f (fe) 
«w 

6l<fc)p('=)+4-a< fc )-a( fc )-p<' i) -(2d-l) 



(fe) 
Po = 

(fe) 
Pi = 

4 fc) =(l-pDa 



(6) 



where we defined — p^ /(l 
malization condition pg' '' + p^ - 



fp( fc )), with the nor- 
4*° = 1. The star 
tionarity condition, summarized in the (^), provides the 
renormalized density vector at each scale, and it couples 
the dynamical properties to the stationary ones. It is 
worth to remark that we do not determine the RG equa- 
tions for p from the coarse graining prescriptions. In fact 
the stationary properties have to be evaluated consider- 
ing the average over many dynamical processes. Thus 
the densities p arc determined from the renormalized dy- 
namical description of the system, namely eq.s @ with 
renormalized parameters. Note also that the RG equa- 
tions are written in the hypothesis of a double time scale 
separation: fire spreading and trees growth are not in- 
teracting. This is expressed by the consistency relation 
p < 1/T(£), where T(£) is the average time needed to 
burn a cluster of characteristic size £. Since T(£) ~ £ 2 
where z is the dynamical exponent, and £ ~ if /p) ~ v we 
have that T(£) ~ (//p^'with v' = vz. Our RG ap- 
proach is then suitable to describe what happens in the 
region of the phase space (p, 9) defined by the condition 
p < 6-"' . 

Given this scheme, we can find the fixed points of the 
renormalization transformation by studying the flow di- 
agram in the phase space of the parameters (p, p, 9) . In 
two dimensions the RG equations show the fixed point 
p* =0,9* =0 and p* = (2/3, 1/3, 0). A complete charac- 
terization of the fixed point is obtained by the linearized 
RG equations around p* ,9* and p* . We find only one 
relevant scaling field that corresponds to an eigenvalue 
greater than 1, which is given by 



A = 



d9^ 



\p*,o 



(7) 



Therefore, we have that the fixed point is repulsive in the 
9 direction, which of course defines the relevant control 
parameter. Since the fixed point is repulsive we can de- 
termine the exponent of the clusters characteristic length 
by using the largest eigenvalue, i.e., v — log 2/ log A. We 
obtain v = 0.73 by using the simple 2x2 cell renormal- 
ization scheme. We can easily improve the results with a 
3x3 cell calculation, which gives v = 0.65, showing that 
the numerical result converges to the right value with re- 
fined renormalization scheme. The exponent v describes 
the divergence of the correlation length by R ~ 8~ v and 
the value obtained is in good agreement with the value 
v ~ 0.58 measured in |qjll]|. In this perspective the pa- 



rameter 



/ 



plays the role of the relevant critical pa- 
rameter as the reduced temperature in the thermal phase 



transitions. For each value of 9, small but finite, the sys- 
tem shows an upper characteristic length in the cluster 
distribution. Only for 9* = the system is critical and 
shows an infinite correlation length. 

The exponent t describing the distribution of fire 
spreading can be obtained as follows. The fires are rep- 
resented by the clusters of connected sites interested by 
a burning process. As in || we define K as the probabil- 
ity that an active relaxation process (i.e. fire) is limited 
between the scales frC 0-1 ) and and it does not extend 
further: 



b< fc > / poo 

A : = / P(r)dr / / P{r)dr 



(8) 



where P(r)dr is the probability to have a burning 
cluster with radius between r and r + dr. In two di- 
mensions with simple scaling arguments we can conclude 
that, if P(s) ~ s~ T and s ~ r 2 (compact clusters), then 
P(r) ~ r 1 ~ 2r . Inserting this expression in eq.(§) we ob- 
tain (in d = 2) r = 1 - (log(l - K))/ (2 log 2). In our 
case K is the probability that at a generic scale (k) all 
the nearest neighbors of a burning tree are empty and 
then K = (1 — p^)*- In the scale invariant regime 
(Pi = pi), K = 0.1975 and then we obtain for r the 
value r ~ 1.16 in excellent agreement with very accurate 
simulations performed by P. Grassberger raj . 

Along the same lines it is possible to compute also 
the dynamical exponent z, and then from scaling laws 
|]llf the remaining critical exponents that characterize 
the FFM model p2[ . Our approach can be naturally ex- 
tended to the one dimensional forest fire. We can follow 
the scheme used previously, and along the same lines it 
is possible to compute the exponent v = 1.0 and r = 1.0 
which recover the exact results for the one dimensional 
case obtained in We summarize our results for the 
one and two dimensional case in Table 1. Details of the 
one dimensional calculation as well as the two dimen- 
sional case will be reported in Ref. fl2]| . 

In conclusion, we introduce a novel RG approach suited 
to study the critical state of the Forest-Fire model. We 
identify the structure of the phase space in which the 
RG transformation is constructed and then we introduce 
the coupling of the renormalization equations with the 
asymptotic stationarity condition of the systems. By 
studying the evolution of the system under scale change, 
we stress the existence of a relevant parameter, 9 7 namely 
corresponding to a repulsive fixed point for the RG equa- 
tions. In this sense, 9 is the critical parameter of the 
model, and the critical state is reached only with a fine 
tuning of 9 to its critical value. The existence of this 
parameter, given by the ratio between the driving rate 
of the system / and the trees growth rate p, places the 
FFM in the field of ordinary critical phenomena. In fact, 
it is worth to stress that in sandpile models the same RG 
analysis || shows that no relevant parameter is present, 
namely the fixed point is completely attractive, and the 
slow-driving of the system does not affect the infinite 
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correlation properties but only the lower cut-off of distri- 
butions. 

It is a pleasure to thank B. Drossel and P.Bak for in- 
teresting discussions. 
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d= 1 


V 


z 


T 


RG Scheme 


1.0 


1.0 


1.0 


Exact results" 1 " 


1.0 


1.0 


1.0 


d = 2 








RG Scheme (2 x 2) 


0.73 


1.17 


1.16 


RG Scheme (3 x 3) 


0.65 


1.02 


1.16 


Simulations* 


0.58 


1.04 


1.15 



In this table we summarize our results for the critical 
exponents compared with exact + p"3| ] or experimental* 
results |0F1. 
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• Fig.l: Configurations of sites at scale M fe ' corre- 
sponding to a green cell at scale b^ k+1 \ 

• Fig. 2: Example of two processes contributing to 
the renormalization of the lightning probability 

• Fig. 3: Example of renormalization of the growth 
parameter . 
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a = 3 



a = 4 



